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Abstract 

Einstein's rate equations are derived from the semiclassical Bloch equations describing the interaction of a classical broadband 
light field with a two-level system. 
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I. INTRODUCTION 

Lorentz's damped and forced harmonic oscillator equa- 
tion pj and Einstein's rate equations |2] are very valuable 
and widely used heuristic models for the description of 
light-matter interaction. Although each of these mod- 
els applies to different situations (roughly speaking, the 
Lorentz equation usually applies when weak or highly de- 
tuned fields are involved whilst Einstein rate equations 
are most adequate for describing resonant interaction of 
atoms with broadband fields), both can be fully justified 
withing the framework of quantum optics in the appro- 
priate limits US HQ. 

Generically speaking, formal derivations from first 
principles of heuristic models are important. This is not 
only for easthetical reasons or completitude arguments 
as these derivations allow to clarify the domain of ap- 
plicability of the heuristic model. The situation of the 
Lorentz and Einstein models is quite different in this 
sense: Whilst formal derivations of the Lorentz equation 
from the semiclassical theory of light-matter interaction, 
i.e. from the optical Bloch equations, can be found in 
several textbooks (see, e.g. UQ), formal derivation of 
the Einstein's rate equations (ERE for short) cannot be 
found in textbooks nor in research articles to the best of 
the authors knowledge This is the most surprinsing 
given the paramount importance of the Einstein model. 

Of course, the above does not mean at all that the con- 
nection between the optical Bloch equations and ERE 
has not been treated: Several quantum optics textbooks 
(see, e.g. discuss ERE and provide derivations 

of the A and B Einstein coefficients. We find it impor- 
tant to remark that whilst the derivation of the A co- 
efficient requires quantization of both matter and elec- 
tromagnetic field (see, e.g. |f|) §], the B coefficient can 
be derived within the standard semiclassical formalism 
(in which only matter is quantized), a derivation that 
provides the same result as that obtained with the fully 
quantized theory p|. Obviously, the values of A and B 
thus obtained verify the relation already given by Ein- 
stein for the ratio A/B. It is also worth noting that the 
semiclassical derivations of the B coefficient make use of 
the assumption of weak field, an approximation that is 
not necessary indeed. 

In this article we derive ERE from the Bloch equations 
with heuristically added spontaneous emission terms as 
we are interested in stimulated processes. The article 
is organized as follows: In section II we present briefly 
ERE; in section III we introduce the optical Bloch equa- 



tions and reduce them to an integro-differential equation 
for the atomic inversion. Then, we derive ERE, first in 
the limit of weak field intensity (section IV) and later 
removing this assumption (section V). We find it con- 
venient to make this distinction because the derivation 
in the weak field limit does not involve certain statisti- 
cal assumptions which are necessary for arbitrary field 
strength. In section VI we discuss the limits of valid- 
ity of ERE and finally, in section VI, we give our main 
conclusions. 



II. EINSTEIN RATE EQUATIONS 

ERE were postulated by Albert Einstein in his fa- 
mous 1917 paper and a clear presentation of them 
can be found, e.g., in Einstein introduced three 

basic processes describing the interaction of a broad- 
band isotropic light field with an atomic two-level sys- 
tem: Stimulated absortion, stimulated emission, and 
spontaneous emission. Einstein assigned to these pro- 
cesses a constant probability per unit time given by 
B 12 W (W21) , B 2 iW (lu 2 i) , and A 2 i, respectively, being 
W (u>2i) the spectral energy density of the light field at 
the frequency u>2i of the atomic transition, and A 2 ±, B\ 2 , 
and B 2 i the Einstein coefficients, which do not depend 
on the field strength nor on time. 

Assuming that all A atoms are in either the upper 
excited state (labelled 2) or in the lower fundamental 
state (1), the time evolution of the population of atoms 
in level 2 (A 2 ) and 1 (Ai = A — N 2 ) are governed by 

A 2 = -A 21 N 2 + W (wax) (B i2 Ni - B 21 N 2 ) , (1) 

Ai - -A 2 , (2) 

assuming that the number of atoms A is large enough for 
the individual absorptions and emissions produce smooth 
temporal changes in the populations (the dot stands for 
time derivative). 

Through the analysis of thermal equilibrium it is eas; 



i it is easy 



to prove that B\ 2 = B 2 \ and A 2 \jB 2 \ = hto 2 
Finally, by introducing the normalized population inver 
sion 

_ A 2 - Ai 



A 



(3) 



(1 > n > — 1) we can rewrite ERE Eqs.flJ and © in the 
simpler form 



dt 



n (t) = -A [n (t) + 1] - 2BW (w 21 ) n (t) , (4) 
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with A = A21 and B = B\ 2 = B 2 i- 
A. Weak field limit of ERE 

For later purpouses it is convenient to consider the 
special case of a weak light field. This can be done by 
writing W (uj 2 i) = eW\ (^21), with e <C 1 and W\ (^21) 
an order one quantity, expanding the inversion in powers 
of e 

00 

n (t) = z m n m (t) , (5) 

m=0 

and writing the equation of evolution for the first orders 
in e. At order e° Eq.0J reads dn /dt — ~A(n a + 1), 
and thus uq = —1 after a transient. Making use of this, 
at order e one gets 

iUi (t) = -A [hi (t) + 1] + 2BW 1 (waO , (6) 
at 

which is the weak field limit of ERE. 



proportional to the A coefficient) since, as it is well 
known, the standard semiclassical theory cannot describe 
spontaneous emission [? ]. 

The optical Bloch equations (|9"|)- (|11[) admit the simpler 
form 

n = -A(n + l)-2i(n*p 21 -n Pl2 ), (13) 

P21 = --Ap^i - iftn, (14) 

where n(r,t) = p 22 (r,t) ~ Pn ( r , t) is the population 
inversion of the atom. 

Now Ea. (|14f> can be integrated formally yielding 

P21 (r, t) = -i / dr ft (r, r) n (r, r) e^^, (15) 
Jo 

where we have taken P21 (r, 0) = for the sake of sim- 
plicity. Substituting this into Eq. l|13fl we get 

n = -A (n + 1) (16) 

- 4 Re / drVl* (r, t) il (r, t) n (r, t) eJ a(t " 4) . 
j 



III. BLOCH EQUATIONS FOR A TWO-LEVEL 
SYSTEM 

Let us consider a light field of the form 

£(r,i)=E(r,t)+E*(r,t), (7) 
E(r,t) =y"d 3 fcE(k)e l(k - r ~ cfct) , (8) 

interacting with a closed two-level atom or molecule, 
whose excited and fundamental states are \2) and |1), 
respectively. The two-level atom has a transition fre- 
quency W21 and electric dipole matrix elements (2 \ p, |1) = 
(1| /i |2) = pz, which have been taken to be real vectors, 
alligned parallel to the cartesian z— axis, without loss of 
generality. In the Dirac picture, and after performing 
the rotating-wave approximation, the semiclassical op- 
tical Bloch equations for an individual atom located at 
point r can be written as 0, 0, 

P22 (r, t) = -Ap 22 - 1 [fl* (r, t) P21 ~ ft (r, t) p 12 ] , (9) 
p n (r, t) = +A P22 + 1 p* (r, t) p 2X - ft (r, t) Pl2 ] , (10) 

hi 0, t) = -~Ap2i - in (r, t) {p 22 - p n ) , (11) 

where p 22 (pn) is the population of the upper (lower) 
atomic level, pi 2 is the slowly varying atomic coherence 
and 

ft(r,i) = ^z-E(r,i)e iW21 *, (12) 

is half the complex Rabi frequency of the light field. The 
effect of spontaneous emission needs to be described phe- 
nomenologicaly (through the damping/pumping terms 



A. Ensemble averaging 

Let us emphasize that Eq. (|16J) is the equation of evo- 
lution for the population inversion for a single atom. As 
we are interested in the average evolution of the whole 
system, the ensemble of atoms, which is the quantity 
whose evolution ERE describe, we next write down the 
evolution equation for the ensemble averaged inversion 

n(t) = (n{r a ,t)), (17) 

where the index a labels each atom and 

1 N 

(/(r Q ,t)) = -^/(r a ,t), (18) 

a=l 

for any arbitrary function / (r a ,t) (N the total number 
ot atoms). The equation of evolution of the ensemble- 
averaged inversion h (t) is obtained from Eq. (|16fl and 
reads 

— n = -A(h + l)-4Re / drC n (t, r) e^-*), 
dt J 

(19) 

C n (t, t) = (ft* (r a , t) ft (r„, r) n (r„, t)> . (20) 

In the next section we derive ERE for the particular case 
of a weak light field, and leave the general case for the 
following section. 

IV. DERIVATION OF THE EINSTEIN'S RATE 
EQUATIONS I: WEAK FIELD 

We start with Ea. (|16|) and follow the same steps as in 
deriving Eq.@ in subsection III Al Assume a weak field 



2 



ft (r, r) = efti (r, r), expand the inversion in powers of 
e, Eq. 10), and find the equation of evolution at order e. 
One easily gets 

hi = -A(n 1 + l)+4Re / drft* (r, t) ft (r, r) e^ A(r - 4) . 

•/ o 

(21) 

Now, performing the ensemble averaging 

— n 1 = -A(n 1 +l) + 4Re / drC (i, r) e^"*), (22) 

Jo 

where C (t, t) is the correlation function 

C(t,r) = (ft*(r a ,t)ft(r ,r)). (23) 

A. Evaluation of the correlation C (t, r) 

In this subsection we calculate the correlation function 
C (t, r) for the special case of an isotropic broadband 
light-field. From Eq. (|2*3^l and definition (fl%)) we have 

1 N 

C(*,r) = -£ft*(r Q ,*)ft(r ,r), (24) 

a=l 

which, under the assumptions of very large N and homo- 
geneous spatial distribution of atoms, can be computed 
as 

C (t, t) = — J d 3 r ft* (r, t) ft (r, r) . (25) 

with V the volume occupied by the sample. Making use 
of Eq. O, Eq. J2HJ) becomes 

C{t ^ = Vtfj d^rEtiv^E^v^e-^-^, (26) 

with E z (r, i) = z • E (r, t). By using Eq. J7J), and assum- 
ing a large enough gas volume as for being possible to 
apply / d 3 r e l ( k '^ k ) r = (2tt) 3 6 (k' - k), which is strictly 
valid only for an infinite volume, Eq. (|26|1 becomes 

C(t,T) = (M^L. e -^i(t-T) J d 3 k |^( k) |2 e icfe( t -r) ) 

(27) 

with E z (k) = z • E (k) 

We now evaluate the integral in Eq. I|27|l by expressing 
the electric vectors as 

E (k) = e, | (k) E ] | (k) + e ± (k) E ± (k) , (28) 

where the polarization unit vectors e|| (k) and ej_ (k) are 
perpendicular to k, and are chosen to be parallel and 
orthogonal, respectively, to the plane defined by k and z. 
With this choice, 



with 9 is the angle defined by k and z (lo|. 

Substituting Eqs. and (gSJl into Eq. |2ZI and ex- 
pressing the integral in spherical coordinates, we obtain 



c(t,T) = ^ 2n \f e -^(t-r) r d r d9sin a e (30) 

Vh 2 J J 

poo 

x dkk 2 \E\\ (k)|V 
Jo 



|2 ick(t- T ) 



We finally assume that |-E(k)| is dependent on k but 
not on the orientation of k, i.e. we assume that the light 
energy at a given k is distributed isotropically in all di- 
rections. In particular, we write 



\E\\ (k)| 2 = \E n {lo)\ 2 , co = ck. 



(31) 



Under this assumption, the angular integrations can be 
done and changing the integration varible from k to uj, 
Eq. (J20J) becomes 



x / dej cj 2 |J5|| (w)| 2 e lw(t - T) . 
Jo 



(32) 



Finally, taking into account that the spectral density of 
energy per unit volume of the light field can be written 
as (see the Appendix) 

W( W ) = ^(27r)V |£||M| 2 , (33) 

we can express the correlation function C (t, r), Eq. I|32|l . 
in terms of W (ui) as 



C(t,T) 



6h 2 e jo 
Introducing [3 = lu — LU21, we get 



e-' iW2l(t - T) / duj W (to) e zuj{t - T) . (34) 



C(t,r) 



6h 2 eo 



d(i W [lo 21 + 13) e lt){t ~ T \ (35) 



where the lower integration limit has been extended from 
—0J21 to —00, which is a very good approximation for op- 
tical frequencies. Hence it turns out that C (i, r) is pro- 
portional to the Fourier transform of the spectral density 
of energy of the light field. 



B. Rate Equations 



Making use of Eq.||33J), Eq. (J221 reads 
d , , , „ , 2u 2 



e|| (k) = sin#, z • ej_ (k) = 0, 



(29) 



x / dp W (wai + 0) K (/?) 

J — oo 

K(p) = Re [ dre-i^-^)^. 



(36) 



(37) 
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Performing the integral K (/3) one gets, 

i _ e -{kA-m)t 



K {(3) = Re ■ 



Re 



(38) 



d(3 W [lo 21 + /3)K^wW (w 2 i) 



which is a Lorentzian of width A centered ar (3 = 0. 
Going back to Eq. (JSSJ), if W (cj 2 i + /?) is a smooth 
function, (specifically, if A [dW (lu 2 i + f3) /d/3]g =0 <C 1) 
we can take W (w 2 i + /3) = W (lu 2 i) and the final result 
reads (f*™ d/3 K (/3) = n) : 

j t ^) = -A(n 1 + l) + ^-W M , (39) 

which compared with Eq. JHJ provides the following value 
for the Einstein B coefficient: 



D 



7T/i 
3ft 2 £Q 



(40) 



which is the correct one for isotropic radiation 



V. DERIVATION OF THE EINSTEIN'S RATE 
EQUATIONS II: ARBITRARY FIELD STRENGTH 

Now we remove the weak light field assumption and 
return to Eq. ((T^I) . Now it is the correlation C„, Eci. (|2*U)l . 
that has to be evaluated. 

We are considering broadband incoherent light fields. 
This implies that the quantities appearing in C n oscil- 
late randomly in time: On the one hand the phases of 
the Fourier components components of the field, E(k), 
are random as the field is incoherent, and on the other 
hand, the values taken by the population inversion of 
a single atom n (r, r) are not correlated with those of 
the fields at this very instant, because the value of the 
inversion depends on previous times. Thus, it seems rea- 
sonable to assume that the random variations of the in- 
version of each individual atom, n(r a ,r), are decorre- 
lated from the random variations of the mutual intensity 
fl* (r a , t) (r a , r). This assumption constitutes our main 
approximation. Thus, we assume 

C n (t, r) » (n* (r Q , t) n (r„, t)) (n (r 0) r)) = C (t, r) n (r) , 

(41) 

with C (t, t) given by Ea.l(55)l. 

Once Ea. (|41|l is assumed, the derivation of ERE is 
easy: substituting Eas. l|41|l and l|35|l into Eq. (|19|) one 

gets 



S fi (*) = -A [n (*) + !] 3h2£o 



Re / drn(r)e-^ (t - T) 
Jo 

(42) 



f -\-oo 



dpW{iu 21 +f3) e Wt - T \ 



which after making the same assumptions as in Sec. IIV Bl 
becomes 



jn (t) = -A [n (t) + 1] - ^-W (u, 21 ) n (t) . 



(43) 



This is the same Eq. (J2J after identifying the B coeffi- 
cient, En.pfijl. 

We see that it is not necessary to assume a weak field 
for deriving the B coefficient with the semiclassical the- 
ory. In spite of that, for arbitrary field intensity it is the 
statistical independence of the flutuations of the inver- 
sion and the mutual intensity that must be invoked. 



VI. ESTIMATION OF THE DOMAIN OF APPLI- 
CABILITY OF THE EINSTEIN'S RATE EQUA- 
TIONS 



The derivation presented above is a bit crude in the 
sense that the approximation carried out in Sect. IIVBI 
(that W (ui) was a smooth enough function of ui) does not 
inform on the domain of validity of the resulting ERE. In 
other words, we have assumed an infinitely broad spec- 
trum for the light field and the natural question is: How 
broad must the spectrum be for ERE be valid? In this 
section we give such an estimation by considering a fi- 
nite spectrum lincwidth. Especifically, we assume that 
W (uj) has a Lorentzian form, which represents many ac- 
tual light sources but can be regarded also as a repre- 
sentation of any bell-shaped distribution. Then we take 



w (w) = w M 



7 2 + (w - luq Y 



(44) 



where loq is the center of the spectrum, which has a width 
7. We note that 



W(wai) = W(u ) 



7 



7 2 + S 2 ' 



(45) 



where 5 — lo 2 i — luq is the detuning between the transition 
frequency and the center of the light field spectrum. With 
this lorentzian form, Eq. 134(1 can be evaluated to yield 



C(t,r) 



7T/1 



7iy(w )e-' i5(t - r) e^ (t ^ T) . (46) 



We now go back to Eq. (|19|) . with approximation 141fl . 
and obtain 



-n (*) = -A [n (t) + 1] - 2 1 BW M (47) 



xRe / dTfi(T)e-(? A+ i +iS P- T \ 
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where we have made use of Eq. (|40|) . We finally evaluate 
the integral by parts repeatedly, yielding 



d , , 



Re 



A [n (t) + 1] - 2jBW {uj ) 
(-if 



(48) 



E 

p=0 



\P+1 



X 71 



V 2 
(P) (*) 



hA + j + i5) 

-ftW (0) e -(^+7+«> 



where the superscript (p) denotes the p— th derivative 
with respect to t. If we now assume that 7 > 4 (the 
spectrum is broad referred to the natural linewidth of 
the atomic transition), the exponentials in Eq. I|48(l can 
be ignored as they are damped out in a time ~ 7 -1 much 
smaller that the characteristic time A -1 of the inversion 
dynamics. Then 



dt 



n (t) = -A [n (t) + 1] - 2-/BW (w 
(-If 



EM*)' 

p=0 



(t)Re 



(7 + i<5) 



P+i ' 



(49) 



(50) 



Clearly, in order that Eq. (|49|l be equivalent to Ein- 
stein equation, Eq. (@J , the terms with p > 1 should be 
negligible. Let us then evaluate the order of magnitude 
of the p— th order, which we denote by O (S p ): 



0(S p )=o(n^ 



Re- 



1 



( 7 + isr 



= O 



(51) 



We thus need to evaluate O (n^). We start with p = 0, 
i.e., (t) = n (t) whose order of magnitude is 1. Thus 
O (n^ )) = 1. We proceed with p = 1, and assume that 
it can be computed by neglecting the S P >\S in Eq. |@SJ. 
We then have 

(t) = -A [n (t) + 1] - 2BT^ (uj 2 i) n (t) , (52) 

where Eq. (|43l has been used. An upper bound 
to the order of magnitude of fS 1 ^ can be estimated 
by using O (h^) = O (ft) = 1, hence O (n^) = 
2A + 2BW (lu 2 i). Similarly we obtain n^> (t) ~ 
- [A + 2iW (w 2 i)] (*), and making use of Eq. ijESjl. 
C(n( 2 )) = [2A + 2.6^(^21 )] 2 , where a factor 2 mul- 
tiplying A has beed added, which amounts to over- 
In general one gets O (fj/ p )) = 
Substituting this expression into Eq. 



estimate 0(n^). 
[2A + 2BW(lu 21 )] p 
B51JI we obtain 



O (S p ) = 7 [2A + 2BW (lu 21 )] p ( 7 2 + S 2 Y 



(53) 



Now we now impose O (S p +i) <C 0(S p ), what directly 
leads to an upper bound to the spectral density, 



2BW (tjgi) 
^7 2 + <5 2 



< 1, 



(54) 



where we have recalled that A <C 7. 

In resume, under condition (|54|l . all 5 p >i can be ne- 
glected and Eq. (|4^|l reduces to the Einstein's rate equa- 
tion Eq. @. Notice that condition (|54|l imposes a con- 
strain on the field spectrum for ERE being applicable: If 
the energy density W (uj 2 {) is very large, the spectrum 
linewidth 7 (or the detuning 8) must be consequently 
large. 



VII. CONCLUSIONS 

We have provided a derivation of the Einstein's rate 
equations from the semiclassical Bloch equation for an 
esemble of two-level atoms or molecules. The derivation 
has been done first for a weak field and we have obtained 
the right expression for the B coefficient. In our deriva- 
tion the necessary ensemble averaging is done explicitly. 
Then, we have generalized the derivation to the case or 
arbitrary field strength by assuming a statistical inde- 
pendence of the flutuations of the atomic inversion and 
the mutual intensity, which is a reasonable assumption 
for an incoherent light field. Finally we have estimated 
which condition the spectrum of such an incoherent field 
must verify for the ERE be valid, Eq. 

Let us finally note that in all the above derivations 
the condition At S> 1 has been invoked. Thus the ERE 
describe the "long term" behavior of the system and their 
solutions are only meaningful in this limit. 



VIII. APPENDIX 

Let us now see how Eq. (|3*2*|l relates to the spectral 
density of energy per unit volume of the light field. The 
energy per unit volume 77 can be computed as 



V = 



£0 
2V 



(ft 



& (v,t) + <?B* (r, t) 



(55) 



By making use of Eq. J7J, writing B(r,t) — B (r, t) + 
B* (r,i) with 

B(r,t)=f d 3 fcB(k)e i(k r ~ cfct) , (56) 

and making use of the well known property (a x b) ■ 
(c x d) = (a ■ c) (b • d) — (a • d) (b ■ c), it can be read- 
ily shown that 



„ = 2^(2tt) 3 / <rk E,ki 



|B(k)| 



Taking into account that 
B(k) = 



k x E (k) 
clkl ' 



Eq. (|57|l simplifies to 

r ? = 2|(2^) 3 | d 3 k\E(k)\ 2 



(57) 



(58) 



(59) 
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Now, making use of Eq. (J2BJ, expressing this integral in |£j| (uj)\ . Finally, from Eq. JSHJ we identify the spectral 
spherical coordinates, and changing the integration vari- density of energy per unit volume as 
able from k to lo = ck, one gets 

" = 8 V? (5hr) 7o ^l^l^l 2 ' ( 6 °) ^M = 8^_(2 7 r) 4 ^|i? || H| 2 . (61) 

where the angular integration has been performed and 
again it has been used the isotropy condition \E± (cu)\ — 



[1] H. A. Lorentz, The Theory of Electrons (Teubner, 
Leipzig, 1909), Chap. 4. 

[2] A. Einstein, Quantum Theory of Radiation, Phys. Z. 18, 
121-128 (1917); an English translation appears in The 
World of the Atom, edited by H. A. Boorse and L. Motz 
(Basic Books, New York, 1966), Vol.2, pp. 888-901. 

[3] P. W. Milonni and J.H. Ebberly, Lasers, (John Wiley & 
Sons, New York, 1988) 

[4] B. W. Shore, The Theory of Atomic Coherent Excitation 
(John Wiley & Sons, New York, 1990). 

[5] R. Loudon, The Quantum Theory of Light (Oxford Uni- 
versity Press, 2000). 

[6] J. N. Dodd, Atoms and Light: Interactions (Plenum 
Press, New York, 1991). 

[7] Let us stress that derivations of rate equations through 
adibatic ellimination of the medium polarization (valid 
in the limit of large decay for the atomic coherence, as 
compared to that of the population inversion, see e.g. 0) 
are not true derivations of Einstein's rate equations. For 
example, in this way one does not obtain the right ex- 
presion for the Einstein's coefficients. For deriving ERE, 



one must take into account the broadband nature of the 
light field and has not to imposse any a priori condition 
on the atomic decay rates. 

[8] This is not completely exact: Although it is true that 
the standard semiclassical theory does not explain spon- 
taneous emission, there exists at least one form of the 
semiclassical theory that can. Self-field Quantum Elec- 
trodynamics, the theory developed by A.O. Barut and 
collaborators during the eighties, is a radiation-matter 
interaction theory that correctly describes spontaneous 
emission without field quantization. See, e.g., [Ij 

[9] A. O. Barut and J. P. Dowling, Self-field quantum elec- 
trodynamics: the two-level atom, Phys. Rev. A 41, 2284- 
2294 (1990). 

[10] Note that the single case in which the referred plane is 
undefined is when k and z are parallel, i.e. — 0, but in 
this case both polarization unit vectors are orthogonal to 
z, in agrement with the result z • e|| (k) = sin# = and 
z ■ e ± (k) = 0. 



6 



